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The cross section for the process e + e — > uiir 7r°7r°7 has been measured in the energy range 
$H ' 1.05-2.00 GeV. The experiment has been performed at the e + e" collider VEPP-2000 with the SND 

Below 1.4 GeV our data are in good agreement with the previous SND and CMD-2 measurements. 
Data on the e + e~ — > cj7t° cross section are well described by the VMD model with two excited p-like 
states. From the measured cross section we have extracted the 7* — > urn transition form factor. It 
has been found that the VDM model cannot describe simultaneously our data and data obtained 
from the ui — > ty° decay. We have also tested CVC hypothesis comparing our results on the 

bjir~ v T decay. 



PACS numbers: 13.66.Bc, 14.40.Bc, 13.40.Gp, 13.25.Gv 



I. INTRODUCTION 



> 

OO ' Experiments at the e+e~ collider VEPP-2000 Q with the SND detector @] started in 2010. The main goals of 
these experiments are a high precision measurement of the total e + e — > hadrons cross section in the energy range up 
1 to 2 GeV and investigation of the vector meson excitations with masses between 1 to 2 GeV/c 2 . In this connection, 
. \ the study of the process 

m 

O ■ e + e~ -> ujtt° -> 7r°7r°7 (1) 

m ; 

is very topical. The process e + e~ — > W7r° is one of the dominant hadronic processes contributing to the total hadronic 
cross section at the center-of-mass (cm.) energy between 1 and 2 GeV. As one of important decay modes of the 
isovector vector states p(1450) and p(1700), it can provide a lot of information about their properties. Moreover, this 
measurement can be used to check the relation between the the e + e _ — > luit cross section and the differential rate 
in the t" — > ujtt~ v t decay following from the conservation of the vector current and isospin symmetry 0] . The SND 
plans to search for electric dipolc decays of the p(1450) and p(1700) mesons to 7r°7r°7 final state, which are important 
for understanding of p(1450) and p(1700) quark structure. The e + e" — > u>ir° is the main background process for this 
search; its precise measurement is the first step in investigating the p(1450) and p(1700) radiative decays. 

In this work the e + e~ — > u>ir° process is studied in the u decay mode to ir°j. Despite the fact that the main lu 
decay mode to 7r + 7r~7r° has probability about an order of magnitude higher, this choice looks reasonable. Unlike the 
47r final state, the W7r° intermediate mechanism is dominant in the 7r°7r°7 final state in the energy range under study. 
This makes it possible to avoid systematic uncertainties due to both the complex procedure of subtracting background 
and taking into account interference between different intermediate mechanisms. 

The process e + e~ — > urn in the ui — > 7r°7 decay mode was first studied in ND experiment Q at VEPP-2M collider. 
The cross section was measured at cm. energies below 1.4 GeV. Later this measurement was repeated by SND Q 
and CMD-2 @ detectors with much higher statistics. Our measurement of the e + e~ — > W7r° — > Tr°n°j cross section 
at VEPP-2000 based on 5 pb" 1 collected in 2010 was published in Ref. 0. 
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The first measurement of the process e + e — > um in the w -> 7r + 7r 7r° decay mode was performed by DM2 0]. 
For a long time this measurement was the only one above 1.4 GeV. Below 1.4 GeV this cross section was measured 
at VEPP-2M by SND 0] and CMD-2 [T3|. 

II. EXPERIMENT 

SND Q is a general purpose non-magnetic detector. Its main part is a spherical three-layer Nal(Tl) calorimeter 
with 560 individual crystals per layer and 90% solid angle coverage. The calorimeter energy resolution for photons 
is <te/E 7 = 4.2%/ yf E 1 {GeV) , the angular resolution ~ 1.5°. There is a tracking system around the collider beam 
pipe based on a nine-layer drift chamber and a one- layer proportional chamber with cathode-strip readout. Outside 
the calorimeter a muon detector consisting of proportional tubes and scintillation counters is placed. An aerogel 
Cherenkov counter located between the drift chamber and the calorimeter is used for particle identification. 

Experiments at VEPP-2000 started in 2010. During 2010-2012 the cm. energy range E = 1.05-2.00 GeV was 
scanned several times with a step of 20-25 McV. The total integrated luminosity collected by SND in this energy 
range is about 40 pb _1 . This work is based on data (27 pb" 1 ) recorded in 2010-2011. 

During the experiment, the beam energy was determined using the measurements of the magnetic field in the collider 
bending magnets. To fix the absolute energy scale, the scan of the ^(1020) resonance was performed and its mass was 
measured. However, possible instability and uncertainties in collider components may lead to a sizable energy bias 
when the beam energy increases from 0.5 to 1 GeV. The uncertainty in the energy setting was investigated in 2012 in 



special runs, in which the beam energy was measured using the Compton backscattcring method [11(. Basing on the 
results of these runs we conservatively estimate the uncertainty in the cm. energy to be 5 MeV. 

III. LUMINOSITY MEASUREMENT 

In this analysis the process 

e + e~ — > 77. (2) 

is used for luminosity measurement. Similar to the process under study |T]) the normalizing process © does not 
contain charged particles in the final state. The selection criteria for the process are chosen in such a manner that 
some uncertainties on the cross section measurement cancel as a result of the normalization. 

For example, in the selection of five-photon events of the process under study we require the absence of a charged 
track in an event. This leads to loss of signal events that contain beam-generated spurious tracks. The probability 
of such a loss may reach several percent and strongly depends on experimental conditions (vacuum pressure in the 
collider beam pipe, beam currents, . . . ). Since the same condition is used for selection of the normalization process, 
the systematic uncertainty associated with beam-generated extra tracks cancels. Moreover, events of both processes 
are selected by the same hardware trigger. Therefore the uncertainty associated with the trigger inefficiency cancels 
too. 

To select events of the process ([2]), the following selection criteria are used: 

• at least two photons and no charged particles are detected, 

• the number of hits in the drift chamber is less than or equal to five, 

• the energies of two most energetic photons in an event are larger than 0.3-E, 

• the azimuth angles of these photons satisfy the condition \\c/>i — <t>2\ — 180°| < 11.5°, 

• the polar angles of these photons satisfy the conditions \6\ + 9 2 - 180°| < 17.2° and 36° < #1,2 < 144°. 

The condition on the number of hits suppresses background from Bhabha events with unreconstructed tracks in the 
drift chambers. We do not expect any significant background from other e + e~ annihilation processes in the energy 
region under study. To estimate possible cosmic-ray background, data recorded during 7.5 hours in a special run 
without beams are analyzed. No e + e~ — > 77 candidates are selected. We estimate that the fraction of cosmic events 
in the sample of selected two-photon events is less than 2 x 10~ 4 and conclude that the cosmic-ray background is 
negligible. 

To calculate the detection efficiency and the cross section of the process ([2]), a Monte-Carlo (MC) event generator 



based on Ref. 12[ is used. The integrated luminosity measured for each energy point is listed in TablcQ] The theoretical 
uncertainty on the cross section calculation is about 1%. The systematic uncertainty on the detection efficiency is 
estimated to be 2%. 
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FIG. 1: The distribution of the 7r°7 invariant mass for selected data events (points with error bars) with E < 1.7 GeV 
(left) and E > 1.7 GeV (right). The curves are the result of the fit described in the text. The dashed line represents the 
linear-background contribution. 

IV. EVENT SELECTION 

At the first stage of the analysis, five-photon events are selected with following criteria: 

• at least five photons and no charged particles are detected, 

• the number of hits in the drift chamber is less than or equal to five, 

• Etot/E > 0.5, where E to t is the total energy deposition in the calorimeter. 

For events passing the preliminary selection, kinematic fits to the e + e~ — > 57 and e + e~ — > 7r°7r°7 hypotheses 
are performed with requirements of energy and momentum conservation and ir° mass constraints for the second 
hypothesis. The goodness of the fits is characterized by the \ 2 parameters: Xs 7 an d X^o 7r o 7 - For events with more 
than five photons, all five-photon combinations are tested and the one with minimal x^-o n oy is used. To select lutt 
candidates the following additional conditions arc applied: 

• Xs 7 < 30 for E < 1.7 GeV and X s 7 < 15 for E > 1.7 GeV, 

• xlo^o 7 - x% < io, 

• at least one of the two 7r°7 invariant masses satisfies the condition |m T o 7 — M u \ < 200 MeV/c 2 , where is 
the cj-meson nominal mass [l3T | . 

The distributions of the 7r°7 invariant mass for 7899 selected data events with E < 1.7 GeV and 331 data events 
with E > 1.7 GeV are shown in Fig. [TJ The ui meson peak is clearly seen in both distributions. Since each event has 
two entries into the histogram, the nonresonant parts of the distributions are determined mainly by signal events. 

V. FITTING THE 7r° 7 MASS SPECTRA 

To determine the number of signal events, the 7r°7 mass spectrum is fitted by a sum of signal (F s ) and background 
(Fb) distributions. The signal distribution is obtained by fitting the mass spectrum for simulated signal events with a 
non-parametric kernel estimation technique fl4j . To take into account a difference between data and MC simulation 
in mass resolution and calibration, the 7r°7 mass for simulated events is smeared and shifted before the fit. The values 
of the smearing Gaussian sigma (<r s ) and mass shift (Am) are deduced from comparison of the w peak width and 
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position in data with the same parameters in simulation. To study energy dependence of a a and Am. the comparison 
is performed in three cm. energy regions below 1.7 GeV. The mass shift and the Gaussian sigma are found to be 
linearly dependent on the energy. The value of Am changes from (1.0 ± 0.5) MeV/c 2 in the energy region 1.05-1.30 
GeV to (3.9 ± 0.5) MeV/s 2 in the energy region 1.5-1.7 GeV. The value of <j s for the same energy regions changes 
from (6.2±1.0) MeV/c 2 to (10.5±0.5) MeV/c 2 . For energies above 1.7 GeV, where statistics are small, a s is obtained 
by a linear extrapolation from lower energies. 

The main sources of background are QED processes such as e + e~ — > 37, 4-f, 5-f, and the hadronic processes 
e + e~ — > ?/7 and e + e~ — > W7r 7r°. Background can arise also from the process e + e _ — > 7r°7r°7 with intermediate states 
other than the W7r°, for example, fa"/. The ir°j mass distribution for the QED, 777 and 7r 7r°7 events is expected to 
be flat and is described by a linear function. The distribution for ujTr°ir° events obtained from MC simulation has 
a complex shape with a wide maximum shifted to the right relative the u peak position. The expected number of 



W7r°7r° events is calculated using experimental data on the e + e~ — > w7t + 7t~ cross section 15j and the isotopic relation 
(r(u;7r + 7r~) = 2cr(w7r°7r°). The background from e + e~ — > ujtt tt is important in the energy range 1.7-1.9 GeV, but 
even here it does not exceed 6% of signal. 

The cross section for the process e + e~ — > p°ir° — > ir°ir°j is estimated from the e + e~ — > pir — > Sir cross section 
measured by BABAR [lfl, and found to be small, below 1 pb. However, because of closeness of the p and u> 
masses, the interference between the W7r° and p°ir° amplitudes can give a sizable contribution to the measured 
e + e~ — > u!tt° — > 7r°7r°7 cross section. The ir°j spectrum for the interference term is peaked at uj mass and practically 
indistinguishable from the spectrum for the lutt intermediate state. Since the phase between ujtt and p°ir° amplitudes 
is unknown, we calculate the maximum possible value of the interference term and use this value as an estimate of the 
systematic uncertainty on the measured lott cross section. The uncertainty due to the interference is estimated to be 
2% at E < 1.55 GeV, then increases to 4.5% at 1.7 GeV, and to 8.0% at 1.8 GeV. Above 1.8 GeV this uncertainty is 
negligible in comparison with the uncertainty of the radiative correction. 

To determine the number of signal events we perform an unbinned extended likelihood ht to the mass spectrum in 
the range \m^o 7 — M^\ < 200 MeV/c 2 . The likelihood function used for E < 1.7 GeV is given as follows: 

N s (l + k s ) 



M 



L = P P (N- N s + N b )P B (M - N; N, k t ) J] F,(ro£o 7 ) 



N s (l + k s ) + N b (l + k b ) 

N h {l + k b ) 
N s (l + k s ) + N b (l + k b ) 



fiK^l .n ■ :\ ■ Et ) (3) 



where N is the number of selected events, M is the total number of entries in the fitted spectrum, N s (N b ) is the 
numbers of signal (background) events, k s (k b ) is the fraction of signal (background) events with two entries per 
event, i.e. events for which masses of both ir°j combinations satisfy the condition |m 7r o 7 — MJ\ < 200 MeV/c 2 , 
kt = {N s k s + N b k b ) I (N s + N b ). The function Pp and Pb arc Poisson and binomial distributions for the total number 
of selected events and the number of selected events with two entries to the fitted spectrum, respectively. The 
parameter k s is calculated using signal MC simulation. It changes from 0.4 at E = 1.1 GeV to 0.2 at E = 2.0 GeV. 
To understand the range of k b variation, we use MC simulation of the e + e~ — > ir°ir°j events at the generator level in 
different models: according to phase space, with the intermediate /o (980)7, /2 (1270)7, an d fo (1370)7 states. In the 
energy range under study the k b value is found to change from 0.1 to 0.5. Therefore, the parameter k b is set to be 
0.3 ±0.2 in the fit (the likelihood function is multiplied by the corresponding Gaussian). The background distribution 
F b is described by a linear function. Above 1.7 GeV a term describing W7r°7r° contribution is added into the likelihood 
function. 

The fit results are shown in Fig.[TJ To obtain the shape of the mass spectrum for signal, the distribution of simulated 
events over energy points is weighted to yield the distribution observed in data. It is seen that our model for signal 
and background describes well the experimental mass spectra. The total numbers of signal and background events 
are 7533 ± 110 and 366 ± 70, respectively, for E < 1.7 GeV, and 282 ± 22 and 49 ± 15 for E > 1.7 GeV. The similar 
fitting procedure is applied in each energy point. The numbers of signal events obtained from the fits are listed in 
Table H 



VI. DETECTION EFFICIENCY AND RADIATIVE CORRECTIONS 



The detection efficiency for the process under study is determined using MC simulation. The simulation includes 
radiative corrections to the Born cross section calculated according to Ref. [T3|- In particular, an extra photon 
emitted by the initial electrons is generated with the angular distribution modeled according to Ref. [lH- The 
detection efficiency e r is evaluated as a function of two parameters: the cm. energy E and the energy of the extra 
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FIG. 2: The E r dependence of the detection efficiency for the process e + e — > cj7r° — > k 7T°7 for E = f.05 GeV (left) and 
E = 2.00 GeV (right). 



photon E r . The E r dependence of the detection efficiency is shown in Fig. [2] for the energy points with minimum and 
maximum energies studied. 



7r°7r°7 is written as 



The visible cross section for the process e + e — > urn 1 

a ms = J £ r (E,xE/2)F(x,E)a(Ey/T^x)d, 



(4) 



where <j(E) is the Born cross section, which one needs to extract from the experiment, F(x, E) is a function de- 
scribing the probability to emit extra photons with the total energy xE/2 [l7T |. Equation (j4]) can be rewritten in the 
conventional form: 



a ms = e(E)a(E)(l + 5(E)), 
where 6(E) is the radiative correction, and e(E) is defined as follows: 

e(E) = e r (E,0). 



(5) 



(6) 



Technically the determination of the Born cross section is performed as follows. With the use of Eq. (j4]) the energy 
dependence of the measured visible cross section is approximated. To do this the Born cross section is parametrized 
by some theoretical model that describes data reasonably well. The fitted model parameters are used to evaluate 
the radiative correction 5(E). Then the experimental values of the Born cross section are obtained using Eq.([5]). To 
estimate the model dependence of the radiative correction, the model parameters are varied in a wide range, with the 
condition that the approximation quality remains acceptable. The fit to the measured cross section is described in 
detail in the Sec. IVIII The obtained values of the radiative correction are listed in Table HI The systematic uncertainty 
associated with the radiative corrections is estimated to be about 1% at E < 1.6 GeV and increases up to 5% at 
E = 1.7 GeV. Above 1.7 GeV the radiative correction becomes large and highly model dependent. For these energy 
points we quote a range of (1 + 5) variation. 

Imperfect simulation of detector response for photons leads to a systematic uncertainty in the detection efficiency. 
To estimate this uncertainty, the data distributions of the most important selection parameters xi 7 and (x^.o 7r o 7 — xi 7 ) 
for signal events are compared with corresponding simulated distributions. The distributions are shown in Fig. [3] for 
events with E < 1.7 GeV. They are normalized to the same number of events. The distributions for background 



events shown in Fig. [3] by the shaded histograms are obtained using data from the 



sideband (66.7 McV/c 2 < 



|m 7r o 7 — M u \ < 200.0 MeV/c ). These background distributions are added to the simulated signal distributions. 
A difference between data and MC simulation is seen in both the distributions. To obtain numerical estimation, 



we change the limits of the conditions on Xs 7 and (x^-o 7r o 1 



xi 7 j 



from 30 to 50 and from 10 to 50, respectively. 



The resulting variation of the measured cross section Jow is found to be (—0.4 ± 0.5)% for the xi 7 condition, and 
(2.9 ± 0.5)% for the (x^o 7r o 7 — Xs 7 ) condition. The obtained values of 5a^ are used to calculate a correction to the 
detection efficiency. The detection efficiency obtained from simulation for our standard selection criteria should be 
decreased by (-2.5 ± 0.7)%. 
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FIG. 3: The xi-y (left) and (x^o^o- — X57) (right) distributions for data events with E < 1.7 GeV (points with error bars). 
The open histogram is a sum of the simulated signal distribution and background distribution. The latter is shown by the 
shaded histogram. The data and the signal+background distributions are normalized to the same number of events. The 
arrows indicate the selection criteria used. 

In SND a photon converted in material before the drift chamber produces a track. Events with converted photons 
are rejected by the selection criteria used. Since the numbers of photons in the final state are different for the signal 
and normalization processes, the data-MC simulation difference in the conversion probability leads to a systematic 
shift in the measured cross section. The conversion probability is measured using — > 777 — > 37 events collected in 
a special run in vicinity of the </>-meson resonance and is found to be (0.97 ± 0.28)% in data and (0.78 ± 0.04)% in 
simulation. The data-MC simulation difference (0.19±0.28)% is consistent with zero. We conclude that the simulation 
reproduces the photon conversion reasonably well. The statistical error (0.3%) of the data-MC simulation difference 
is taken as an estimate of the systematic uncertainty on the photon conversion. The corresponding uncertainty on 
the measured cross section is estimated as 3 x 0.3 = 0.9% (the contributions from two of five photons cancel due to 
normalization) . 

As it was discussed earlier, a some part of the data events contains beam-generated spurious tracks and photons. The 
effect of extra tracks cancels due to normalization to two-photon events. The presence of extra photons also changes 
the detection efficiency, but differently for the signal and normalization processes. To take into account this effect 
in MC simulation, beam-background events recorded during experiment with a special random trigger are merged 
with simulated events. The detection efficiencies obtained using simulation with and without merged background 
are compared. It is found that the presence of extra photons do not influence the number of selected normalization 
events. The detection efficiency for the signal process increases by 0.3-1.3% depending on experimental conditions. 
Unfortunately, the random-trigger events were not recorded on a regular basis during the 2010-2011 experiments. 
Therefore, we conservatively estimate the correction due to extra photons to the signal detection efficiency to be 0.8% 
with a systematic uncertainty of 0.5%. 

The corrected values of the detection efficiency e(E) are listed in Table [I] The statistical error on the detection 
efficiency is negligible. A systematic uncertainty from the sources discussed above is 1.2%. A nonmonotonic behavior 
of e{E) as a function of the cm. energy is due to variations of experimental conditions, in particular, due to change 
of the number of dead calorimeter channels during data taking period. 

VII. RESULTS AND DISCUSSION 

The Born cross section for e + e~ — > uir° — > 7r°7r°7 obtained using Eq.([5]) is shown in Fig. [4] in comparison with 
results of previous measurements. The numerical values are listed in Table |TJ The quoted errors on the cross section 
are statistical and systematic. The sources of systematic uncertainty are summarized in Table [TTJ The total systematic 
uncertainty is 3.4% in the energy range E < 1.55 GeV and 4.5% in the energy range 1.55 < E < 1.6. Above 1.6 GeV 
the uncertainty increases due to the model dependence of the radiative correction. 

Our data are in good agreement with the measurements H, i] performed by SND and CMD-2 at the VEPP- 
2M collider at energies below 1.4 GeV, but significantly (by 20-30%) exceed the DM2 data @|. The DM2 data 
were obtained in the 7r + 7r~7r 7r° mode and have been rescaled using the ratio of the w — > ir°j and U) — > Sir decay 
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TABLE I: The cm. energy (E), integrated luminosity (IL), detection efficiency (e), number of selected signal events (N s ), 
radiative-correction factor (1 + 5), measured Born cross section (a). For the cross section the first error is statistical, the second 
is systematic. 



E, GeV 


IL, nb- 1 


e, % 


N s 


1 + 8 


a, nb 


1.050 


358 


35.5 


104 ± 11 


0.903 


90 ± 10 ± 03 


1.075 


545 


36.5 


176 ± 15 


0.913 


0.97 ± 0.08 ± 0.03 


1.100 


845 


36.0 


297 ± 17 


0.921 


1 06 ± 06 ± 04 


1.125 


518 


35.9 


220 ± 16 


0.928 


1.28 ± 0.09 ± 0.04 


1.150 


412 


37.8 


178 ± 13 


0.934 


1.23 ± 0.09 ± 0.04 


1.175 


539 


36.8 


231 ± 17 


0.939 


1.24 ± 0.09 ± 0.04 


1.200 


1058 


36.6 


489 ± 24 


0.943 


1.34 ± 0.06 ± 0.05 


1.225 


550 


37.8 


265 ± 19 


0.947 


1.35 ± 0.10 ± 0.05 


1.250 


435 


37.9 


187 ± 17 


0.950 


1.19 ± 0.11 ± 0.04 


1.275 


495 


37.0 


254 ± 21 


0.953 


1.46 ± 0.12 ± 0.05 


1.300 


1278 


37.6 


673 ± 35 


0.956 


1.47 ± 0.08 ± 0.05 


1.325 


522 


38.2 


279 ± 24 


0.959 


1.46 ± 0.12 ± 0.05 


1.350 


554 


38.1 


292 ± 24 


0.962 


1.44 ± 0.12 ± 0.05 


1.375 


574 


38.1 


302 ± 24 


0.966 


1.43 ± 0.11 ± 0.05 


1.400 


1012 


37.9 


578 ± 33 


0.970 


1.55 ± 0.09 ± 0.05 


1.425 


598 


38.1 


363 ± 26 


0.977 


1.63 ±0.12 ±0.05 


1.450 


427 


38.3 


221 ± 19 


0.985 


1.37 ±0.12 ±0.05 


1.475 


599 


38.4 


291 ± 21 


0.995 


1.27 ±0.09 ±0.04 


1.500 


1939 


39.0 


996 ± 40 


1.007 


1.31 ±0.05 ±0.04 


1.525 


487 


38.2 


245 ± 19 


1.021 


1.29 ±0.10 ±0.05 


1.550 


543 


38.2 


228 ± 16 


1.038 


1.06 ±0.08 ±0.04 


1.575 


505 


37.9 


170 ±17 


1.063 


0.83 ±0.08 ±0.04 


1.600 


814 


38.1 


232 ± 19 


1.100 


0.68 ±0.06 ±0.03 


1.625 


505 


37.9 


139 ± 15 


1.161 


0.63 ± 0.07lo.o4 


1.650 


473 


36.9 


96 ± 10 


1.262 


0.44 ± 0.04±° .03 


1.675 


454 


37.0 


74 ±11 


1.429 


0.31 ± 0.05lo'o2 


1.700 


698 


30.3 


70 ± 10 


1.704 


0.19±0.03l^i 


1.725 


502 


30.6 


22 ±6 


2.0-2.3 


0.06 ± 0.02±g;gg5 


1.750 


503 


29.2 


25 ±6 


2.4-3.3 


0.06 ± O.OllHi 


1.775 


521 


28.7 


22 ±6 


2.8-5.0 


0.033 ± 0.008l^o? 


1.800 


727 


27.9 


33 ±7 


3.4-9.0 


0.023 ± 0.005i^o 5 


1.825 


477 


28.2 


7±3 


4-15 


0.004 ± 0.002l^ 08 


1.850 


400 


26.7 




5-24 


0.002 ± 0.002l^o 05 


1.870 


631 


26.2 


19 ±6 


5-27 


0.005 ± 0.002l°;° 16 


1.890 


577 


26.3 


24 ± 5 


6-31 


0.006 ± O.OOllHo? 


1.900 
1.925 
1.950 
1.975 


553 
555 
406 
460 


24.9 
24.6 
23.3 
24.0 


I2±t 

utl 

It? 


6-40 
5-68 
5-63 
5-39 


0.004 ± 0.002l°;°" 

0.006 ± 0.002jl ( o : o^ 
n nni +0.002+0.001 
u.uui_ .001-0 .001 

0.009 ± 0.0041^07 


2.000 


536 


23.3 


5±l 


4-39 


0.006 ± 0.0031^05 



TABLE II: The systematic uncertainties on the measured cross section from different sources and the correction to the detection 
efficiency. The total uncertainty is the sum of all the contributions in quadrature. 



Source 



Systematic uncertainty, % 



Correction, % 



Luminosity 

Selection criteria 

Photon conversion 

Beam background 

Radiative correction (E < 1.6 GeV) 

Interference with p°7r° (E < 1.6 GeV) 

Total 



2.2 
0.7 
0.9 
0.5 
1 

2-3.6 
3.4-4.5 



-2.5 
0.8 
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C 

c 1.6 
o 

1.4 

8 1.2 
o 

1 

0.8 
0.6 
0.4 
0.2 




• SND2013 
A SND 2000 

★ CMD-2 
□ DM2 








E (GeV) 



FIG. 4: The cross section for e + e 



7r°7r°7 



measured in this work (circles), and in SND [f| (triangles), CMD-2 



(stars), and DM2 [8j] (squares) experiments. Only statistical errors are shown. The curve is the result of the fit to SND 2000 
and SND 20f3 data described in the text. 



probabilities jl3| . 

The cross section measured in this work is fitted together with the SND data obtained in experiments at VEPP- 
2M The Born cross section is described by the following formula Q 



a(E) = 



47TCT 

~E^~ 



.9 pLO-K 

~17 



D 



D, 



D 



A^" 



p" 



D 



P f {E)B(uj -> 7r° 7 ), 



(7) 



where a is the fine structure constant, g puJ ^ is the p — > lot: coupling constant, f p is the 7* — > p coupling constant 



calculated from the p 



decay width, D Pi (E) 



E — iEY Pi {E), m Pi and T Pi (E) are the mass and 



width of the resonance pi, B(lo — > is the w — > 7r°7 branching fraction. The first term in Eq. ([7]) describes the 
contribution of the p(770) resonance, the second and third represent the /9(1450) and p(1700) contributions. The 
fourth term is added to study a model dependence of the measured cross section due to possible existence of a broad 
p-likc resonance with mass m p <» > 2 GeV /c 2 or due to nonresonant contribution into the Born cross section near 2 
GeV. The parameters Ai are the ratios of the coupling constants Ai = g p 4tJ7r / 'g ' pul7V ■ f p /f Pi , <fii and </>2 are the phases 
of the p(1450) and p(1700) amplitudes relative to the p(770) amplitude. For p'" contribution, the phase is assumed to 
be equal to zero. The energy dependence of the phase-space factor Pf{E) is calculated using the MC event generator 



for signal events. For an infinitely narrow u> resonance, Pf{E) 
energy dependence of the p(770) width is described as follows: 



1/3 • q^j, where q u is the uj meson momentum. The 



T P (E) - T p (m p ) 



E 



qn{m p ) 



q 2 
LZtt 



(8) 



where q w (E) = -J (E/2) 2 — m 2 , is the ir~ mass. For p(1450) and p(1700), the energy-independent widths are 
used. 

The data are fit with free parameters g pul7r , A±, A2, A3, M p i, M p », <f>i, and fa- The values of the p' and p" widths 
are fixed at PDG values [TJ|: V p , = 400 MeV, I> = 250 MeV. The parameters M p „, and V p „, are set to 2.3 GeV 
and 400 MeV, respectively. The fit result is shown in Fig. [4] The fitted parameters are used to calculate the values 
of the radiative corrections listed in Table HI To estimate the model dependence of the radiative correction, the fit is 
performed with M p > fixed at different values from the range 1.4-1.6 GeV and with A3 either free or fixed at zero. The 
quality -P(x 2 ; v ) of these fits, where v is the number of degrees of freedom, varies from 5 to 20%. 
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TABLE III: Fit parameters obtained. 



Parameter 


Model 1 


Model 2 


o GeV -1 


15 6 ± 3 

1 'J ■ W — 1 — V / . > J 


17.4 ± 0.1 


A x 


0.26 ± 0.01 


0.11 ± 0.001 


A 2 


0.060 ± 0.006 


= 


M p , , MeV 


1491 ± 19 


1488 ± 3 


I>, MeV 


= 400 


321 ± 4 


M p », MeV 


1708 ± 41 




I>, MeV 


= 250 




0i, deg. 


168 ± 3 


121 ± 2 


^2, deg. 


10 ± 7 




X> 


56.8 / 52 


118.6 / 54 



To study the contributions of the p(1450) and p(1700) resonances, we restrict the energy range to E < 1.9 GeV. 
This reduces the model uncertainty due to a possible nonresonant contribution or the p'" resonance. The parameter 
A 3 is set to zero. The fit results are presented in Table Hill Two models have been studied, with non-zero and zero p" 
contribution. In the fit with the second model the parameter F p i is left free. Such a fit with only one p excitation was 
performed in Rcf. Q and gave a reasonable description of CMD-2 and DM2 data. It is seen that our more precise 
data cannot be described by the model with one excited p state. The fitted value of the parameter Ai is used to 
calculate the products of the branching fractions 



B( Pi 



B{ Pl 



e+e 



) = 



o>, {m Pi )m 2 p . 
12tt 



(9) 



where 



4ira 2 



~17 



Ai — — 



Pf(E) 



(10) 



is the cross section for the process e+e — > p- L — > um without interference with other p-like resonances. The results 
are following 



B(p' -> e+e-) ■ B(p' -> wtt°) = (5.3 ± 0.4) x 10 



-6 



B{p" -> e+e-) ■ B{p" -> wtt°) = (1.7 ± 0.4) x 10~ 6 . (11) 

It should be noted that there is no generally accepted approach at the moment for describing the tail of the p(770) 
resonance above 1 GeV and shapes of broad resonances like p 1 and p" . The excitation curves of the three resonances 
p, p' and p" overlap; their amplitudes strongly interfere with each other. As a result, a small change of the resonance 
shape can lead to significant shifts in fitted resonance parameters. So, the results obtained with our very simple model 
using energy-independent p' and p" widths can be considered only as rough estimates of the resonance parameters. 

The cross section for e+e~ — s- lutt can be expressed in term of the 7* — > ujit transition form factor F UJ7V1 {q 2 ) (T9I [20I] . 
where q is the four-momentum of the virtual photon: 

Attc^ 2 

°^o(E) = —\F^(E 2 )\ 2 P f {E). (12) 

This form factor is measured also in the cj — > ir Q e+e~ [2l], [2^] and u> — > n°p+p~ [2^, [24[ decays. The value of the 
form factor at q 2 = is related to the u> — > ir°j partial width: 

r(c^ A) = f^ 7 3 IW0)l 2 , (is) 

where P 1 is the decay photon momentum. Using Eqs. ([7]), (fT2"|) . (fT3"]) , and data from Table HTT1 we calculate T(oj — > 
7r°7) for the model 1 to be 0.88 ± 0.05 MeV. For such a simple model the agreement with the experimental value 
r(w 7T°7) = 0.703 ± 0.024 MeV looks reasonable. 

Figure O shows the normalized transition form factor squared {\F U T T1 {q 2 ) / ^^(O)! 2 ) measured in this work and in 
Ref. |5 together with most precise data from omega decays obtained in NA60 experiment (2~H . The curve represent 
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FIG. 5: The 7* — ► u)n° transition form factor. The points with error bars represent data from this work (circles), Ref. 0] 
(triangles), and Ref. [24| (squares). Only statistical errors are shown. The curve represent the result of model prediction with 
the parameters listed in Table Hill for the model 1. The dashed curve shows the p(770) contribution. 



the results of model prediction with the parameters listed in Table Hill for the model 1. The dashed curve shows the 
p(770) contribution. We conclude that it is hard to describe data from e + e~ annihilation and the to — > 7r°/i + /i - decay 
simultaneously with our model based on vector meson dominance (VMD). 

The conserved vector current (CVC) hypothesis establishes a relation between the charged hadronic current in the 
r decay and the isovector part of the electromagnetic current. So, the e + e _ — > ujir° cross section can be related with 
the spectral function of the r~ — > lq-k~v t decay (V^) 3 : 

= —^-Vum(E). (14) 

The comparison of the e + e~ — > uiir — > 7r°7r°7 cross section measured by SND with the cross section calculated under 
the CVC hypothesis from the r~ — > ujtt^v t spectral function measured in CLEO experiment [25[ is presented in 
Fig. |B1 It is seen that the e + e~ and r data are in reasonable agreement. 

A more quantitative test of the CVC hypothesis can be made by comparing the measured r — > untv T branching 
fraction with the value calculated from the e + e~ — > ujn a cross section according the formula [1, HH 

r(r" ^- Vt ) = ^J^ 3 J <Z 3 K - + 2q 2 )cr U)n0 ( q )d q , (15) 

where \V u d\ is the Cabibbo-Kobayashi-Maskawa matrix element, m T is the r lepton mass, Gf is the Fermi constant. 
We integrate the fitted curve shown in Fig. 2] and obtain the value of the product r(r~ — > um~ v T )B(u> — > 7r°7) = 
(3.68 ±0.04 ±0.13) x 10~ 6 eV. Using the values of the r lifetime and B(uj — > tt ^) we calculate the branching fraction 

B(t~ um~v T ) = (1.96 ± 0.02 ± 0.10) x 10" 2 . (16) 

which is in good agreement with experimental value (1.95 ± 0.08) • 10~ 2 obtained as a difference of the PDG [l3[ 
values for B(t~ — > ujh~v T ) and B{t~ — > ujK~v t ). We conclude that the CVC hypothesis for the W7r system works 
well within the reached experimental accuracy of about 5%. 
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FIG. 6: The cross section for e + e~ — > ivir — > 7r°7r°7 measured in this work (circles) and in Ref. 0] (triangles). Only statistical 
errors are shown. The cross-section data shown by squares was calculated under CVC hypothesis from the spectral function of 
the r — » luhv t decay measured in CLEO experiment (25|. The curve is the result of the fit to SND 2013 and SND 2000 data. 



VIII. SUMMARY 



The cross section for the e + e~ — > urn — > 7T 7T 7 process has been measured with the SND detector at the VEPP- 
2000 e + e _ collider in the energy range of 1.05-2.00 GeV. This is the most accurate measurement of the e + e~ — > unr° 
cross section between 1.4 and 2.0 GeV. Below 1.4 GeV our data agree with the earlier measurements of the same 
reaction performed at the VEPP-2M collider with the SND Q and CMD-2 Q detectors. A significant disagreement 
is observed with DM2 data Q in the energy range 1.3-2.0 GeV. 

Data on the e + e~ —> ljtt° cross section are well described by the VMD model with the three p-like state: p(770), p' 
and p". However, the full data set on the 7* — > W7r° transition form factor including data from both e + e~ annihilation 
and oj decays, in particular, oj — > n°p + p~ [24| . cannot be described by such a simple model. 

We have also tested CVC hypothesis comparing the energy dependence of the e + e~ — > urn cross section with the 
spectral function for the t~ — > ojtt~i> t decay, and calculating the branching fraction for this decay from e + e~ data. 
We have concluded that the CVC hypothesis for the lutt system works well within the experimental accuracy of 5%. 
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